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The measure of event-by-event fluctuations in high-energy heavy-ion collisions corresponds to 
the second moment of the fluctuating quantity distribution of interest. It is shown that the measure 
based on the third moment preserves the properties of $ but those related to the higher moments 
do not. In particular, only the second and third moment measures are intensive as thermodynamic 
quantities. The $2— and $3— measure of p±— fluctuations are computed for the hadron gas in 
equilibrium and the results are analyzed in context of the experimental data. 
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Large acceptance detectors allow one for a detailed analysis of individual collisions of heavy-ions at high-energies. 
Due to hundreds or even thousands of particles produced in these collisions, variety of statistical methods can be 
04 ' applied. There are several interesting proposals Jl]-^] to use the fluctuation measurements as a potential source of 
information on the collision dynamics. However, one faces a problem how to disentangle the 'dynamical' fluctuations 
from the 'trivial' geometrical ones due to the impact parameter variation. The latter fluctuations are large and 
dominate the fluctuations of all extensive event characteristics such as multiplicity or transverse energy. Using the 
' fluctuation (or correlation) measure <f>, which was introduced in our paper Jl[, resolves the problem in a specific 
way. By construction, <E> is exactly the same for nucleon-nucleon (N-N) and nucelus-nucelus (A-A) collisions if the 
A-A collision is a simple superposition of N-N interactions. Consequently, <!> is independent of the centrality of A-A 
collision in such a case. On the other hand, <E> equals zero when the inter-particle correlations are entirely absent. 
& ' The $— measure can be applied to the fluctuations of kinematical quantities such as the event energy or transverse 
1 , momentum and to the fluctuations of event chemical composition |(|[7j . 

The NA49 Collaboration plans to study the chemical fluctuations in a near future S while the data on the transverse 
momentum fluctuations have been already published PJTc|]. The value of $ Pi in the central Pb-Pb collisions at 158 
GeV per nucleon has appeared to be smaller than expected. It has been also claimed jl0| that the correlations, which 
^ , are of short range in the momentum space, are responsible for the nonzero positive value of $ p± being observed. 
The result has been widely discussed In particular, our calculations of & p± in the equilibrium hadron gas 

show |L4| that the positive value of $ Pi appears due to the boson statistics of pions. When the hadronic system at 



freeze-out is identified with the pion gas, the calculated slightly overestimates the experimental value llOj] but, 
as discussed here, the inclusion of the pions which come from the resonance decays removes the discrepency. 

The $— measure corresponds to the second moment of the fluctuating quantity, say the event transverse momentum. 
Recently, it has been suggested |l7j to use the higher moments in an analogous way. However, the authors of |l7| ] 
have not realized that the fluctuation measures based on the higher moments, except that of the third one, do 
not posses a key property of $ which has been mentioned above. Namely, &nn — &aa if the A-A collision is a 
simple superposition of N-N interactions. When treated as thermodynamic quantities, the second and third moment 
measures are intensive while the higher moment ones are not. The aim of this note is to substantiate the comment 
and to discuss usefulness of the third moment measure. We focus our attention on the p±— fluctuations which have 
been already studied experimentally |j|,[l0| . 

Let us first introduce the measure. One defines the single-particle variable z = x — x with the overline denoting 
averaging over a single particle inclusive distribution. The event variable Z, which is a multiparticle analog of z, is 
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defined as Z A = J2iLi( x i ~ x )i where the summation runs over particles from a given event. By construction, (Z) = 0, 
where (...) represents averaging over events. Finally, the measure is defined in the following way 



dcf 



(N) 



(1) 



There is an obvious generalization of the definition (|l|) suggested in |17|. Namely, 



(N) 



(2) 



The fact that <& n = 0, when no inter-particle correlations are present, is evident Jl],[l7|]- We are now going to show 
that $3 as $2, m contrast to <!>„ with n > 3, possesses another nontrivial property which is so useful in the data 
analysis. Namely, $2 and $3 are independent of the source number distribution if the particles originate from several 
identical sources. Then, $2 and $3 are independent of the impact parameter if the A-A collision is a superposition 
of N-N interactions. Let us prove this property. 

-Pi (A) is the normalized distribution of X = Y]^-y %i, when the particles come from the single source. If we have 
k sources distributed according to pk, the A— distribution reads 

p(x) = zZp k j dXl dXk p i( Xfc ) 5 ^ x - ( Xl + + Xfe )) • 



k=l 



The moments of P(X) are 



<A™) 



dcf 



dxx n p(x) = (-i) n — HQ) 



(3) 



where the generating function T equals 

/°° i k 

dXe^ x P{X) = Y,Pk [H(Q) 
k=l 

with 

H{Q) = J dX e iQX Pi(X) . 

Using eq. (0), one computes the first five moments of P(X) as 
(A) = (k) (X)i , 

(A 2 ) = (fc)(A 2 ) 1 + (fc(fc-l))(X)?, 

(A 3 ) = (k) (X 3 )i + 3{k(k - 1))(A 2 ) 1 (A) 1 + (k(k - l)(fc - 2))<A) 3 , 

(X A ) = (fc)(A 4 )i + 4(fc(fc - 1))(A 3 ) 1 (A) 1 + 3<fc(fc - 1))<A 2 } 2 

+3(fc(fc - l)(fc - 2))(A 2 ) 1 (A) 2 + (k(k - l)(fc - 2)(fc - 3))(X)i , 

(A 5 ) = (k){X 5 ) 1 + 5(k(k - 1))(A 4 ) 1 (A) 1 + W(k(k - 1))(A 3 ) 1 (A 2 ) 1 
+7(k(k - l)(k - 2))(A 3 ) 1 (A) 2 + 9(k(k - l)(k - 2))(X 2 ) 2 i(X)i 
+7(k(k - l)(k - 2)(k - 3)}(A 2 )!(A) 3 + (k{k - l)(k - 2){k - 3)(fc - 4)}<A)f , 



where 



/oo 
dX X n Pi (A) and (k n ) d = ^ k 

k=i 



Pk 
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Applying these formulas to the variable Z and taking into account that by definition (Z) = (Z)\ = 0, we get 



(Z 2 ) = (k)(Z 2 } 1 , 
(Z 3 ) = (k){Z 3 ) 1 , 

{z*) = (k){z i ) 1 +mk-i)){z 2 )i, 

(Z 5 ) = (k)(Z 5 ) 1 + 10(fc(fc - l))(Z 3 ) 1 (Z 2 ) 1 . 

Since (N) = (k){N)-L, one finds that 

(Z 2 ) {Z 2 ), (Z 3 ) {Z 3 ), 



(N) (NU ' (AO (NU 

but analogous formulas do not hold for (Z 4 ) and (Z 5 ). Instead, 

(Z*) (Z*)! {k(k-l)){Z 2 )l 
(N) {N)i + (k) <7V)i ' 

{Z 5 ) (Z 5 ), , (fctfc-i))^ 3 )^ 2 ) 
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(N) (N)i (k) (N)i ' 

Therefore, $2 and $3 are independent of the source number distribution while $4, $5 and, obviously, 4> n with n > 5 
do depend onp^. The inclusive distribution, which determines z n , is, of course, independent of the source distribution. 
The above results also show that $2 and $3 are intensive quantities, i.e. they are independent of the system size, 
while with n > 3 are not. Indeed, when the source number is fixed, (k l ) = k l and one observes that only $2 and 
$3 do not dependent on k. Let us note here that the independence of k is, in principle, a weaker requirement than 
the independence of pk- 

The $2 — measure is sensitive to the fluctuations or correlations of various origin. For example, it acquires a finite 
value, which is positive for bosons and negative for fermions, due to the quantum statistics |l4[ . The correlation 
between the particle multiplicity and their kinematical characteristics also influences $2 0]- The energy-momentum 
conservation and presence of the collective motion introduces additional inter-particle correlations. Thus, one con- 
cludes that the nonvanishing value of $2 signals the existence of the correlations in the system but it does not explain 
their origin. In such a situation, $3 seems to be very useful. Indeed, simultaneous measurements of $2 and $3 might 
help to identify the fluctuations which dominate in the system. For this purpose one should theoretically estimate 
contributions of various correlations to $2 and $3 

In our paper we have discussed how to compute $2 in the ideal quantum gas. Now, we are going to extend 
these calculations to the case of $3 . For comparison, we also present here the earlier published results on $2 ■ At 
first, the energy fluctuations arc considered. Therefore, the single particle variable x is identified with the particle 
energy E. Then, one immediately finds that 

< 4 » 

where the single particle average energy is 



1 f d 3 p E 



while the particle density p equals 



p J (2tt) 3 \- x eJ iE ± 1 

d 3 p 1 
(^j 3 \- 1 ef 3E ±l ' 



(5) 



f3 = T^ 1 is the inverse temperature; A = denotes the fugacity and p, the chemical potential; E = y/ m 2 + p 2 with 
m being the particle mass and p its momentum; the upper sign is for fermions while the lower one for bosons. 
Since Z = U — NE, where U is the system energy, (Z 2 ) and (Z 3 ) are computed as 



3 



(Z 2 ) = 5 



(Z 3 ) 



d 2 9 2 ^ 



c^ 2 



2£A 



<9/?<9A 



S(V,T,A) , 



(6) 



d 3 



3E w x §\ +iE2 ~k{ x ik 



S(V,T,A) , 



where H(V", T, A) is the grand canonical partition function jjJI defined as 



N a 



with V denoting the system volume and the index a numerating the system quantum statesQ As well known 
the grand canonical partition function of the quantum ideal gas equals 



\nZ(V,T,\) = ±gV J ^ In [l ± A e^ E ] 



(7) 



with g being the number of the particle internal degrees of freedom. After a rather lengthy calculation, one finds 



(Z 2 ) 



d 3 p 



and 



(N) pj (2tt) 3 
(Z 3 ) 1 f d 3 p 



(E - Ef 



A"V B 



(N) P J (2tt) 



(E-E) 



3 A^e^fA^e^Tl) 
[\- 1 eP E ±l) 3 



(8) 



(9) 



As expected, $2 and $3, which are given by the formulas (|4|,|3||9|) , are intensive thermodynamic quantities, i.e. they 
are independent of the system volume. We also note that $2 and $3 are independent of g. One observes that the 
sign of $2 is definite i.e. $2 < for fermions, $2 > for bosons and $2 = in the classical limit (A -1 3> 1) The 
sign of $3 is not definite but $3 still vanishes for the classical particles. 

When the particles are massless and their chemical potential vanish (A = 1), the calculations can be performed 
analytically to the end. Then, eqs. (H,@,0) give 



' ! 0.40/ 



where the upper case is for fermions and the lower one for bosons. For m = p = the bosonic $3 appears to be 
(logarithmically) divergent due to the singular character of the function (eP E — 1)~ 3 at E — > 0. 

One immediately modifies eqs. (f|,||,^) for the case of the transverse momentum. The respective equations read: 



d 3 p 



(PJ 



1 



X -l e 0E ± 1 > 



(10) 



(Z 2 



(N) 



d 3 p 
(2tt) 3 



Ox -P±_f 



{\- 1 ef }E ±l) 2 1 



(11) 



(Z 3 ) 



d 3 p 



Ox -P±) 



3 A-V^A-^Tl) 



(N) pj (2tt) 3 v ^ ^ j (\- 1 ef 3E ±l) 3 
where p± = psinO with p = |p| and O being the angle between p and the beam (z) axis, and 



Px = ~ 
P 



d 3 p p± 
J2nf X- x e^ E ± 1 



1 The formulas from [jwj analogous to Q and (ffl) are erroneously written but the final results are correct. 



(12) 
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In Figs. 1 and 2 we present with dashed lines the $2 — and $3— measure of p± — fluctuations in the ideal pion gas. 
The pions are, of course, massive (m, = 140 MeV), so $2 and $3 are found numerically from eqs. ( [io| , [TT| , p^| ) . The 
calculations are performed for several values of the pion chemical potential. In the chemical equilibrium p — 0. As 
seen, $2 is positive but $3 is negative. At T = 200 MeV and p = 70 MeV, $3 experiences a rapid growth. This 
happens because the first term from eq. (||) changes the sign from positive to negative at T = 200 MeV. 

It is a far going idealization to treat a fireball at freeze-out as an ideal gas of pions. A substantial fraction of the 
final state pions come from the hadron resonances. These pions do not 'feel' the Bose-Einstein statistics at freeze-out 
and consequently the values of $2 and $3 should be significantly reduced. We estimate the role of resonances in the 
following way. The spectrum of pions, which originate from the resonance decays, is not dramatically different than 
that given by the equilibrium distribution JlSj ] . Therefore, we treat the fireball at freeze-out as a mixture of 'quantum' 
pions - those called 'direct' - and the 'classical' pions which come from the resonance decays. Since the weighting 
functions in eqs. dlfl,|ll|,0) are all equal to \e^ l3E in the classical limit, the formulas analogous to (Efl,OJl^) are 



_ 1 

v n — 



d 3 p 
(2^)3 



(Pa 



1 



\~ l eP B - 1 



+ A r e 



-/3E 



(13) 



(Z 2 



(N) 



d 3 P 
(2^)3 



— \2 



(P± - P± 



(A -1 e^- B - l) 2 



+ A r e 



-pE 



(14) 



(N) 



with 



d 3 p 
(2^)3 



(pj_ - p ± 



d 3 p 
J^nj 3 



A~ 1 e^(A- 1 e^ + 1) 
(X-i e P E - l) 3 



P± 



1 



\~ 1 eP E - 1 



(15) 



P : 



d 3 p 
(2^)3 



\~ l eP E - 1 



+ \ r e^ E 



The parameter A r is chosen is such a way that the number of 'classical' pions equals to the number of pions from 
the resonance decays. Thus, A r is temperature dependent. In the actual calculations, we have taken into account 
the lightest resonances: p(770) and a->(782) which give the dominant contribution. The life time of p, which is 1.3 
fm/c, is not much longer than the time of the fireball decoupling and some pions from the p decays can still 'feel' the 
effect of Bose statitistics. Therefore, the contribution of p to the 'classical' pions is presumably overestimated in our 
calculations. Since we neglect the heavier resonances and weakly decaying particles, which also contribute to the final 
state pions, the two effects partially compensate each other. In any case, our calculations show that the resonances 
do not change the values of <I>2 and $3 dramatically in the domain of temperatures of interest. 

In Figs. 1 and 2 the solid lines represents $2 — and $3— measure which include the resonances. The chemical 
potentials of p and 10 are assumed to be equal to that of pions. As seen, the role of the resonaces is negligible at 
the temperatures below 100 MeV but above this temperature the resonances reduce the fluctuations noticeably. As 
already mentioned, $2 — measure of p±— fluctuations has been experimentally measured in the central Pb-Pb collisions 
by the NA49 collaboration. The first result has been published as $2 = 0.7 ± 0.5 MeV || but the value of $2 is 
increased to 4.6 ±1.5 MeV when the two-track resolution effect is properly taken into account |l0| . If we identify the 
system freeze-out temperature with the slope parameter deduced from the pion transverse momentum distribution 
T = 180 MeV J2(J. Then, the value of $2, which is read out from Fig. 1 for p = 0, equals 15 MeV for no resonances 
and $2 = 8.7 MeV when the resonances are included. The temperature is significantly reduced if the transverse 
hydrodynamic expansion is taken into account. The freeze-out temperature obtained by means of the simultaneous 
analysis of the single particle spectra and the two-particle correlations is about 120 MeV ^0|. Then, the value of 
$2 for p — equals 6.5 MeV for the case of no resonances and $2 = 5.6 MeV when the resonances are included. 
The latter number agrees perfectly well with the mentioned above experimental value. This strongly supports the 
claim jl0| that the short range correlations due to the Bose-Einstein statistics of pions play a dominant role in the 
hadronic system produced in central heavy-ion collisions. However, it would be very interesting to check whether the 
experiment also confirms our prediction on $3 which is presented in Fig. 2. As seen, $3 = —12.3 MeV for T = 120 
MeV and p — when the resonances are taken into account. 
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Let us close this paper with a technical remark. When the measure is applied to the real data or simulated 
events, it is rather inconvenient to use the formula (Q) because then one has to process the data twice; in the first 
run one evaluates the inclusive average x and then computes the moments of Z and z. To avoid the double data 
processing one can use the formula derived in Jl2[ which is 

_/(x 2 ) 2(x)(xn) (xy(Ni) \ 1/2 f(x 2 ) (xy\ 1/2 

2 \(N) (AO 2 (AO 3 J \(N) (N)*J ' 1 ) 

where the event variable X2 is defined as X2 d = J2iLi x 1- The expression of $3, which is analogous to (|l6|), reads 



f{X 3 ) 3(X)(X 2 N) 3{X) 2 (XN 2 ) (X) 3 {N 3 ) \ 1/3 f (X3) 3(X 2 )(X) 2(X) 3 \ 1/3 
V (N) (N) 2 + (TV) 3 (N) 4 J V (AO (AO 2 + (N) 3 ) 



withx 3 d ^Eti4 

We conclude our study as follows. The $3— measure, which is based on the third moment of the fluctuating quantity 
distribution, preserves the advantageous properties of $2 while the higher moment measures do not. Simultaneous 
usage of $2 and $3 may help in identifying the origin of correlations observed in the final state of heavy-ion collisions 
at high-energies. In particular, the measurement of $3 of p±— fluctuations can decisively confirm that the dominant 
correlations in the central collisions are those of the quantum statistics. 



I am very grateful to Marek Gazdzicki for fruitful discussions and critical reading of the manuscript. 
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Figure Captions 



Fig. 1. $2— measure of p±— fluctuations in the hadron gas as a function of temperature for four values of the chemical 
potential. The resonances are either neglected (dashed lines) or taken into account (solid lines). The most upper 
dashed and solid lines correspond to \i = 70 MeV, the lower ones to \i = 0, etc. 

Fig. 2. $ 3 — measure of p±— fluctuations in the hadron gas as a function of temperature for four values of the chemical 
potential. The resonances are either neglected (dashed lines) or taken into account (solid lines). The most upper 
dashed and solid lines correspond to fi — 70 MeV, the lower ones to fj, = 0, etc. 
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